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Abstract 

Background: Areal bone mineral density is predictive for fracture risk. Microstructural bone parameters evaluated at the 
appendicular skeleton by high-resolution peripheral quantitative computed tomography (HR-pQCT) display differences 
between healthy patients and fracture patients. With the simple geometry of the cortex at the distal tibial diaphysis, a 
cortical index of the tibia combining material and mechanical properties correlated highly with bone strength ex vivo. The 
trabecular bone score derived from the scan of the lumbar spine by dual-energy X-ray absorptiometry (DXA) correlated ex 
vivo with the micro architectural parameters. It is unknown if these microstructural correlations could be made in healthy 
premenopausal women. 

Methods: Randomly selected women between 20-40 years of age were examined by DXA and HR-pQCT at the standard 
regions of interest and at customized sub regions to focus on cortical and trabecular parameters of strength separately. For 
cortical strength, at the distal tibia the volumetric cortical index was calculated directly from HR-pQCT and the areal cortical 
index was derived from the DXA scan using a Canny threshold-based tool. For trabecular strength, the trabecular bone 
score was calculated based on the DXA scan of the lumbar spine and was compared with the corresponding parameters 
derived from the HR-pQCT measurements at radius and tibia. 

Results: Seventy-two healthy women were included (average age 33.8 years, average BMI 23.2 kg/m 2 ). The areal cortical 
index correlated highly with the volumetric cortical index at the distal tibia (R = 0.798). The trabecular bone score 
correlated moderately with the microstructural parameters of the trabecular bone. 

Conclusion: This study in randomly selected premenopausal women demonstrated that microstructural parameters of the 
bone evaluated by HR-pQCT correlated with the DXA derived parameters of skeletal regions containing predominantly 
cortical or cancellous bone. Whether these indexes are suitable for better predictions of the fracture risk deserves further 
investigation. 
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Introduction 

In addition to low bone mass, osteoporosis was characterized in 
the WHO Technical Report of 1994 [1] by the micro architectural 
deterioration of bone tissue. Predominantly driven by sex steroid 
hormone levels changing with aging, the structures of the long 
bones are altered by enhanced endosteal resorption and are only 
in part compensated by periosteal apposition [2]. This change 
translates into altered micro architecture of the trabecular and 
cortical bone in vivo as assessed by high resolution imaging 
techniques, such as high-resolution peripheral quantitative com- 
puted tomography (HR-pQCT) [3] . For the cortex, trabeculariza- 



tion of the cortical bone with aging has been suggested [4]. 
Accordingly, ex-vivo micrographs of the femoral cortex demon- 
strated increasing porosity and thinning of the cortex by 
coalescence of the intracortical pores [5] . 

Different studies have demonstrated the direct relationship 
between the mechanical properties of specific anatomical regions 
and their architecture ex vivo [6,7] and in silico [8-10]. Even if the 
differences in the areal bone mineral density (aBMD) measure- 
ments by dual-energy X-ray absorptiometry (DXA) at the 
standardized sites were inexistent or mild, these structural 
parameters could distinguish between individuals with or without 
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prevalent low velocity fractures [11,12]. The standard regions of 
interest (ROI) used in the clinically applied HR-pQCT measure- 
ments at the distal radius and tibia were formerly focused on the 
trabecular bone [13], but meanwhile the quantification of the 
cortical bone derived indices, including the cortical porosity 
(Ct.Po) based on automated segmentation from HR-pQCT 
images, has been implemented [14]. 

A low aBMD is predictive of the individual fracture risk in 
postmenopausal women when measured at the hip, lumbar spine, 
or radius [15,16]. The best predictive ability was observed when 
the measurements were performed at the same skeletal sites that 
sustained a fracture later. Whereas the tibial epiphysis (T-EPI) 
consists predominandy of trabecular bone, the tibial diaphysis (T- 
DIA) is predominantly composed of cortical bone. The authors 
previously showed that, in a cohort of randomly selected 
postmenopausal women, aBMD measured by DXA at both tibial 
subregions was predictive for non-vertebral fracture risk [17]. T- 
scores were generally higher at the T-DIA than at the T-EPI 
suggesting that the T-DIA might be suitable for identifying other 
contributors of an individual's fracture risk, including parameters 
related to the shape and size. 

Because of the growing share of cortex from distal to proximal 
the ratio of trabecular and cortical bone is highly variable at the 
tibial metaphysis [18]. A more proximal region of interest with 
thicker and more homogeneous cortices, such as the T-DIA, has a 
simple geometric structure, similar to a hollow cylinder. The 
combination of the aBMD at the T-DIA and the structural 
parameters derived from DXA scans and confirmed by HR- 
pQCT can improve the prediction of local bone strength ex vivo 
[19]. It is unknown whether this combination of material and 
structural competence can be confirmed by both of the techniques 
in a standardized procedure in vivo. 

The trabecular bone score (TBS) is a DXA-derived parameter 
that focuses on the cancellous architecture. This score quantifies 
the gray-level texture out of the DXA scan of the lumbar spine 
[20]. TBS was shown not to correlate or to correlate weakly with 
the aBMD in postmenopausal women ranging from 0.2 to 0.58 
[21-23]. In epidemiological studies, TBS was significantly lower in 
the women with incidental fractures compared to the controls 
[24]. Even improvements in the TBS could be observed in 
postmenopausal women under antiresorptive treatment in a 
randomized controlled trial [22]. While the correlation of this 
texture parameter to the parameters of microarchitecture, e.g., the 
degree of connectivity, the trabecular number and separation, was 
performed ex vivo by micro-computed tomography [25], little is 
known if these correlations can be confirmed by the HR-pQCT 
measurement of these parameters at the appendicular skeleton in 
vivo. 

The aim of this study was to evaluate whether and to what 
extent the DXA-derived cortical and trabecular structural 
parameters correlate with the bone microstructural parameters 
assessed using HR-pQCT in premenopausal women. 

Materials and Methods 

Study participants 

In December 2012, study invitations were mailed to a total of 
485 women selected from commercially available address listings 
categorized by sex, age, and the region of Davos, Switzerland. 
Mobile, non-institutionalized, not pregnant (confirmed by a 
negative pregnancy test) randomly selected women between 20- 
40 years old were eligible. The exclusion criteria were a known 
bone disease or immobilization (relative or absolute immobiliza- 
tion was defined as unable to stand or walk at least one hour per 



day). In total, 72 women were included in this study, sixty-nine of 
these women were white. The left upper and lower extremities 
were examined or the nonfractured contralateral extremity in 
cases of prior fracture on the left side (Fig. 1). All of the 
measurements occurred at the AO Research Institute Davos, 
Switzerland. 

Ethical Statement 

The study was approved by the local ethics committee for 
medical sciences (Ethics Committee of the Canton of Zurich, 
Switzerland, registration number 2012-0464). All of the partici- 
pants provided prior written informed consent and were mone- 
tarily compensated for their efforts. 

DXA Measurements 

DXA densitometer measurements (Hologic Discovery C™, 
Hologic, Bedford, MA, USA) were performed with a mobile 
device provided by the Department of Osteoporosis, University 
Hospital Berne, Switzerland. aBMD was measured at the lumbar 
spine (LI -4; LS), the femoral neck (FN) and total hip (TH), the 
distal radius (Rad 1/3), and at the two customized subregions of 
the distal tibia (T-EPI, T-DIA). The detailed standardized 
procedures for the tibial BMD measurements by DXA were 
previously published [26]. The ROI is defined as the 120 mm 
height and 129 mm width area that starts 10 mm above the top of 
the ankle joint space. T-EPI corresponds to the distal 40 mm of 
the ROI, and T-DIA corresponds to the proximal 40 mm of the 
ROI (Fig. 2). aBMD is expressed in grams of hydroxyapatite (HA) 
per square centimeter of the projected area. The obtained aBMD 
values of the femoral neck were expressed as the standard 
deviation (SD) of the mean value of the young healthy population 
(T-scores) derived from the National Health and Nutrition 
Examination Survey in the United States of America (NHANES 
III) [27]. The scans were performed according to the manufac- 
turer's guidelines, and the analyses were conducted according to 
the recommendations of the International Society of Clinical 
Densitometry. Spine phantom scanning was routinely performed 
on a daily basis for quality control. A least significant change at the 
95% confidence level of 0.020 g/cm 2 for the aBMD at the total 
hip was achieved as the reproducibility index for the technician 
involved. 

DXA-Derived Parameters of Cortical and Trabecular Bone 

Assuming a hollow cylindrical shape of the T-DIA, the DXA- 
derived cortical thickness (Ct.Th) was calculated using the 
DICOM images of the DXA scan at the distal tibia. DXA gray- 
scale images were imported into Matlab™ programming envi- 
ronment (MATLAB™, MathWorks Inc., Natick, MA, USA). The 
region of interest was automatically extracted for all data sets. 
Image data was segmented using an edge detection filter based on 
the algorithm proposed by Canny [28]. Canny threshold was kept 
constant at 0.05 resulting in a binary image showing the inner and 
outer borders of the cortex. A Butterworth filter was applied for 
smoothing the edges (Fig. 3). Cortex diameter (DIA) was defined as 
mean pixel distance between the outermost edges. Ct.Th was 
calculated as mean distance between outer and inner edges of the 
medial cortex projection, because the individual variation of Ct.Th 
is wider at the medial cortex than at the lateral cortex, i.e., the 
cortex closer to the fibula (own observation, unpublished data). 
Values were transformed to millimeter by isometric scaling (2.0 
pixels/mm). Estimating the bony cross section as ideal tube the 
polar moment of inertia (pMOI) was computed according to 
formula I. 
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Invitations sent 
485 



Did not respond 
382 



Responded 
103 



Did not meet inclusion criteria 
31 



Included 

72 



LS: 
Hip: 
Tibia: 
Radius: 



DXA 

72 

L71 / R 1 
L 70 / R 2 
L 66 / R 6 



HR-pQCT 

Tibia 22.5 mm: L69/R3 
Tibia 60 mm: L69/R3 

Radius 9.5 mm: L58/R6 



Figure 1. Consort flow diagram. 

doi:1 0.1 371 /journal.pone.0088946.g001 




Figure 2. Regions of interest at the left distal tibia by DXA (left) and the volumes of interest as 3D reconstructions by HR-pQCT 
(right). 

doi:1 0.1 371 /journal.pone.0088946.g002 



PLOS ONE | www.plosone.org 



3 



February 2014 | Volume 9 | Issue 2 | e88946 



Bone Parameters by HR-pQCT & DXA in Healthy Women 



Figure 3. ROI of a representative DXA data set of the tibial 
diaphysis. The cortex edges as derived from image filtering are 
depicted as blue (outer diameter) and red lines (intramedullary 
channel). 

doi:1 0.1 371 /journal.pone.0088946.g003 



pMOI = 



DIA 

IT' 



-Ct.Th 



A new parameter was introduced as cortical index for diaphyseal 
bone strength (CI), denned as product of a material (aBMD) and a 
structural component (pMOI). Because this CI is calculated from 
the two-dimensional DXA data, it is called the areal CI (aCI) in 
formula II. 



aCI = aBMD* pMOI 



(II) 



TBS was evaluated in the identical ROI as those used for the LS 
BMD. The dimensionless TBS values were calculated as the mean 



of the individual measurements for vertebrae L1-L4 from the 
DXA using TBS iNsight® Software (version 1.8.2; Med-Imaps, 
Bordeaux, France). 

HR-pQCT Measurements 

HR-pQCT (XtremeCT, Scanco Medical AG, Bruttisellen, 
Switzerland) measurements were performed by immobilizing the 
left forearm and the left distal tibia in a carbon fiber shell and 
scanning as previously described [29-31]. In addition to the 
standard ROI at the radius (radius 9.5 mm) and at the tibia (tibia 
22.5 mm), a third customized ROI was scanned starting at the 
distal tibia 60 mm above the endplate (tibia 60 mm) to increase 
the amount of cortex measured (Fig. 2). A stack of 1 10 parallel CT 
slices was acquired at the distal starting points of all of the sites in 
the proximal direction using 60 kVp and a 126 mm field of view 
reconstructed across an image matrix size of 1536x1536, with a 
voxel size of 82 |J.m. These provided a 3D image of 9.02 mm in 
the axial direction. The attenuation data were converted to the 
equivalent hydroxyapatite (HA) densities. The manufacturer's 
phantom was scanned daily for quality control. 

The bone morphological parameters were evaluated using the 
standard protocols provided by the manufacturer. The parameters 
of interest were the total volumetric bone mineral density (vBMD), 
the bone volume fraction (BV/TV), the trabecular bone density 
(Tb. BMD), the trabecular thickness (Tb.Th), the trabecular 
number (Tb.N), the trabecular separation (Tb.Sp), the inhomoge- 
neity of the trabecular network, defined as Tb.l/N.SD, and 
connectivity density of the trabecular network (Conn.D), the 
cortical bone density (Ct.BMD) and the cortical thickness (Ct.Th). 
The cortical porosity (Ct.Po) and the polar moment of inertia 
(pMOI), as a structural index of resistance to torsion, were 
calculated using an automated image processing chain [14,32]. As 
an analogue to the calculation from DXA, the vBMD at both ROI 
of the tibia multiplied by the pMOI leads to the volumetric cortical 
index (vCI) as a potential parameter of cortical strength. Scans 
with severe movement artifacts, such as a blurry trabecular area or 
skewed or discontinued cortical faces, were excluded. Based on the 
categories (1 =best, 5 = worst), seven forearm scans of grade 3 
and higher were excluded [33]. Additionally, a single forearm 
measurement by HR-pQCT was lost for technical reasons. 



Table 1. Baseline characteristics. 



n = 72 


mean 


SD 


Age (years) 


33.8 


5.5 


Height (cm) 


167.6 


6.2 


Weight (kg) 


65.2 


12.1 


BMI (kg/m2) 


23.2 


4.1 


Calcium intake (mg per day) 


816.0 


402 


n % 


Nulliparas 


42 


58.3 


Smoker total 


19 


26.4 


former 


10 


13.9 


current 


9 


12.5 


Family history for osteoporosis 


22 


30.6 


BMI < 18.5 kg/m 2 


3 


4.2 


Ever glucocorticoid use 


6 


8.3 



doi:1 0.1 371 /joumal.pone.0088946.t001 



Table 2. Areal BMD and the structural 


derivates by DXA. 




DXA 




mean 


SD 


aBMD 


FN 


g/cm 2 


0.853 


0.119 


TH 


g/cm 2 


0.952 


0.116 


LS 


g/cm 2 


1.039 


0.138 


Rad 1/3 


g/cm 2 


0.712 


0.046 


T-DIA 


g/cm 2 


1.362 


0.108 


T-EPI 


g/cm 2 


0.750 


0.096 


Structural derivates 


Ct.Th (T-DIA) 


mm 


3.8 


0.58 


pMOI (T-DIA) 


mm 4 


1.9298 


0.64 


aCI 


g*cm 2 


2.62 


0.856 


TBS 




1.399 


0.087 


doi:1 0.1 371 /joumal.pone.0088946.t002 
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Table 3. Volumetric BMD and the structural parameters by HR-pQCT. 





HR-pQCT 




radius 
9.5 mm 




tibia 22.5 mm 




tibia 60 mm 






Unit 


mean 


SD 


mean 


SD 


mean 


SD 


total vBMD 


mgHA/cm 3 


307.6 


53.3 


303.2 


50.1 


550.3 


79.3 


BV/TV 


% 


13.7 


2.8 


14.1 


2.9 


7.6 


2.6 


Cortical parameters 


Ct.BMD 


mgHA/cm 3 


846.2 


44.7 


874.8 


33.2 


997.9 


20.9 


Ct.Th 


mm 


0.760 


0.151 


1.229 


0.221 


2.521 


0.267 


pMOl 


mm 4 


7694 


2228 


41104 


11140 


20031 


4967 


CtPo 


% 


2.8 


1.0 


4.3 


1.2 


1.2 


0.5 


vCI 


gHA*cm 


0.234 


0.072 


1.231 


0.352 


1.079 


0.222 


Trabecular parameters 


Tb.BMD 


mgHA/cm 3 


164.4 


34.2 


169.7 


34.5 


91.2 


31.2 


Tb.N 


1/mm 


1.944 


0.246 


1.912 


0.312 


1.476 


0.424 


Tb.Th 


mm 


0.070 


0.009 


0.074 


0.011 


0.052 


0.011 


Tb.Sp 


mm 


0.453 


0.073 


0.464 


0.094 


0.690 


0.247 


Inhomogeneity 


mm 


0.187 


0.038 


0.209 


0.053 


0.373 


0.177 


Conn.D 


mrrT 3 


3.300 


0.816 


3.39 


1.01 


1.020 


0.570 



doi:1 0.1 371 /joumal.pone.0088946.t003 



Statistical Analysis 

The analyses were conducted with PASW statistical software 
(version 21.0; SPSS™, Chicago, IL). The two-sided p values 
<0.05 were considered to indicate statistical significance. The 
descriptive data were presented as the mean ± standard deviation 
(SD). After checking the variables of interest for normality 
(Kolmogorov-Smirnov test), the correlations between the struc- 



tural parameters were expressed as R Pearson's coefficient. To 
compare correlations, Fisher's transformation was applied. 

Results 

Baseline Characteristics 

Seventy-two healthy premenopausal women were included in 
the study. The anthropometric baseline characteristics are shown 



5 l 



R = 0.798 



0 -I 1 . 1 1 1 

0.5 0.75 1 1.25 1.5 1.75 

Volumetric strength 
(vCI in gHA*cm) 

Figure 4. Correlation between the aCI (DXA) and the vCI (HR-pQCT) at the distal tibia. 

doi:1 0.1 371 /journal.pone.0088946.g004 
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in Table 1. The mean age was 33.8 years; the mean height and 
BMI were 167.6 cm and 23.2 kg/m 2 , respectively. Forty-two 
women were nulliparas, and 9 of 72 women were current smokers. 
The mean calcium intake was 816 mg per day. Thirty-eight 
females had prevalent fractures, mainly in the forearms (18), 
phalanges (14) and ankle (6). Twenty-two women or 30.6 % 
reported a potentially positive family history for osteoporosis, i.e. 
any self-reported major osteoporotic fracture (clinical spine, 
forearm, hip or shoulder fracture) among first-/second-grade 
relatives. As shown in Table 2, the femoral neck BMD was 
0.853 g/cm 2 , corresponding to a mean NHANES T-score of 0.04 
(± 1.08) SD. All of the areal and volumetric density parameters at 
all of the ROI were normally distributed (Tables 2 and 3). While 
the magnitude of extent of the vBMD was similar at the standard 
ROI of the radius 9.5 mm and the tibia 22.5 mm, the one at the 
tibia 60 mm was higher because of the higher portion of the 
cortical compartment. The Ct.Po was lower at the tibia 60 mm 
compared to the tibia 22.5 mm (1.2±0.5 % vs. 4.3± 1.2 %). 

Correlation Analysis Between the Cortical Parameters 

For the cortex at the distal tibia, a significant correlation (R = 
0.74, p<0.01, Table 4) was found between the CI measured at the 
tibia 22.5 mm by DXA and HR-pQCT, indicating that the bone 
strengths assessed by both techniques were predictive of each 
other. The correlation coefficient was even higher at the tibia 
60 mm (Fig. 4), where the portion of cortical bone is more 
substantial, but without reaching statistical significance (p = 0.3). 

Correlation Analysis Between the DXA-derived and HR- 
pQCT Parameters of the Trabecular Bone 

The TBS of LI -4 correlated significandy (p<0.01) with the 
following parameters at the radius 9.5 mm and the tibia 22.5 mm: 
positively with the Tb.BMD, BV/TV, Tb.N, Conn.D, and 
negatively with the TbSp, inhomogeneity (Table 4). Correlations 
between the TBS and the Tb.Th, the total vBMD, the Ct.BMD, 
the Ct.Th or the Ct.Po were absent or inconsistent between the 
different skeletal sites. Among these healthy women, the TBS 
correlated with the aBMD at the lumbar spine and femoral neck 
(R = 0.660 and 0.532). 

Discussion 

The CI and TBS represent DXA derived parameters of the 
cortical and cancellous bone. This study showed for the first time 
that the corresponding structural parameters of the bone evaluated 
by HR-pQCT correlate with these DXA derived parameters in a 
cohort of randomly selected premenopausal women. In clinical 
practice, the aBMD is the only bone-specific parameter currendy 
used to assess the fracture risk. By integrating the cortical and 
cancellous bone of a defined skeletal area, DXA scans provide 
bone mass of the projected area as expressed as the aBMD. If 
specific skeletal sites revealed a substantial amount of cortical or 
cancellous bone, these structural parameters could be quantified. 
The distal tibial diaphysis offers simple geometry similar to a 
hollow cylinder and a high percentage of cortical bone. Ex vivo 
material and structural characteristics of the bone at the T-DIA 
were previously shown to be predictive of the bone strength 
assessed by biomechanical testing [19]. As shown in the present 
study, the determination of the pMOI on a DXA scan of the tibia 
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This study has some limitations. All of the measurements were 
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necessarily translate into findings of other cohorts, e.g., postmen- 
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included the number of randomly selected healthy women and the 
consistency of the study results of the hip aBMD to the data 
obtained by the NHANES III, which demonstrates the balance of 
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In this cohort of randomly selected healthy women, micro- 
structural parameters of bone assessed by the 3D technique were 
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scans. This approach may add information to improve the 
prediction of individual fracture risks in a clinical setting. 
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